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In order to improve thermoacoustic efficiency, the thermoacoustic devices have been developed from 
standing wave devices to traveling wave devices. Actually, the acoustic field in practical thermoacoustic 
devices is neither a pure traveling wave nor a pure standing wave. The thermoacoustic effect is the hybrid 
effect of traveling wave component and standing wave component. Therefore, the thermoacoustic effect 
of traveling-standing wave will be study in this paper. Firstly, the thermoacoustic conversion perfor¬ 
mance of the traveling-standing wave are analyzes qualitatively by combining the thermoacoustic con¬ 
version performance of the traveling wave with those of the standing wave. Then, based on the basic 
thermoacoustic formulas, the influence of the parameters of the acoustic field and the regenerator’s 
structure on the thermoacoustic conversion is analyzed, and the optimum condition for the thermoacou¬ 
stic conversion is discussed. The results are consistent with the qualitative analysis. Additionally, our the¬ 
oretical results also show a good agreement with the experimental data [Biwa et al. Phys Rev E 
2004;69(6):066304(6)], which indicates the validity of the analysis in this paper. Furthermore, the anal¬ 
ysis in this paper further provides a more intensive understanding of these experimental results. The con¬ 
clusions obtained in this paper are significant to guide for the design of new thermoacoustic devices. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

Thermoacoustic researches have made great progress over dec¬ 
ades, and benefited greatly from the publication of thermoacoustic 
conversion theory. In the early time, various standing wave devices 
[1,2] were built based on thermoacoustic theory. In the standing 
wave device, gas parcels with cp z = (2 n + 1) x 90° realize the conver¬ 
sion between heat and acoustic energy through irreversible thermal 
contact between the working gas and the stuff in the regenerator, 
where cp z is the leading phase of oscillating pressure pi to oscillating 
velocity iq. However, the purposely imperfect thermal contact on 
which the standing wave thermoacoustic conversion relies limits 
its efficiency. With the development of thermoacoustic study, the 
traveling wave thermoacoustic device was proposed by Ceperley 
[3] in 1979. In such engine, gas parcels work at the leading phase 
(p z = 2n x 90° to realize the reversible thermoacoustic conversion 
in theory. However, the reversible thermoacoustic conversion real¬ 
ized by reversible thermal contact requires infinitesimal hydraulic 
radius of the regenerator. This will lead to infinite viscous dissipa¬ 
tion which will weaken the thermoacoustic gain and efficiency 
strongly. So the finite hydraulic radius is adopted in the real travel¬ 
ing wave device, which makes the oscillating temperature of gas 
parcels lag the oscillating displacement of gas parcels. It indicates 
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that the phase difference between the oscillating temperature and 
the oscillating pressure deviates from (2n + l) x 90° because of 
the thermo-relaxation. Under such condition, for improving ther¬ 
moacoustic gain and efficiency, the traveling-standing wave phase 
is necessary to match the phase deviation due to the thermo-relax¬ 
ation. Biwa et al. [4] demonstrated thermoacoustic energy conver¬ 
sions by making full use of the acoustic field induced in the 
resonator. He claimed that, in order to achieve high gain and effi¬ 
ciency, one has to choose an optimum cp z which can make both 
the traveling wave component and the standing wave component 
contribute to the amplification of acoustic intensity. However, the 
leading phase cp z can be chosen in the region of -180° ^cp z ^ 180° 
for the traveling-standing wave, then how to choose the optimum 
(p z for the regenerator becomes a new problem. So far, most 
researchers still focus their particular attentions on the traveling 
wave [5,6] or the standing wave [7-10], and few published litera¬ 
tures have studied on the traveling-standing wave field. 

Actually, the acoustic field in practical thermoacoustic devices 
is neither a pure traveling wave nor a pure standing wave. The 
thermoacoustic effect is the hybrid effect of the traveling wave 
component (TWC) and the standing wave component (SWC). The 
functions of the SWC are energy storage and energy conversion, 
while the functions of the TWC are energy transfer and energy con¬ 
version. In the thermoacoustic engine (TAE), the acoustic energy 
produced by the SWC can be transferred by the TWC, and the 
acoustic energy produced by the TWC also can be deposited with 
the SWC state. In the thermoacoustic cooler (TAC), the acoustic 
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energy consumed by the SWC can be supplied by the TWC, and the 
acoustic energy consumed by the TWC can also be compensated 
from the stored acoustic energy in the SWC. 

So, the thermoacoustic effect of traveling-standing wave will be 
discussed in this paper. The thermoacoustic conversion perfor¬ 
mance of the traveling-standing wave will be analyzed by combin¬ 
ing the thermoacoustic conversion performance of the traveling 
wave with those of the standing wave. Then, based on the basic 
thermoacoustic formulas, the influence of the parameters of the 
acoustic field and the regenerator’s structure on the thermoacou¬ 
stic conversion is analyzed, and the optimum condition for the 
thermoacoustic conversion is discussed. Additionally, some 
numerical results comparing with the experimental data [4] will 
be provided, which will validate the analysis of this paper. 

2. Qualitative analysis 

The thermoacoustic effect is the hybrid effect of the TWC and 
SWC in the practical thermoacoustic system. So the conversion 
capabilities of the traveling wave and the standing wave will be 
analyzed firstly. And then, the thermoacoustic conversion perfor¬ 
mance of the traveling-standing wave by combining the thermoa¬ 
coustic conversion performance of the TWC with those of the SWC 
will be discussed. 

2A. Standing wave mode 

As shown in Fig. 1, the regenerator is located on the left side of 
the pressure antinode (PAN). Taking the rightward direction as the 
positive direction of x, the leading phase cp z of the acoustic pressure 
Pi relative to the velocity iq equals to -90° in the regenerator. In x 
direction, the heat power Q 2 transferred by the regenerator can be 
expressed as [1] 

0-2 =-■^ns k T m pp- i Ui(T-1), ( 1 ) 

and the produced acoustic power W 2 is expressed as [1] 

w 2 = In5 /( Ax^^( Pl ) 2 (r-i) > (2) 

where n is the perimeter, Ax is the length of the regenerator, co is 
the angular frequency, p is the thermal expansion coefficient, Sk and 
b v are the thermal penetration depth and the viscous penetration 
depth respectively, c p is the isobaric specific heat, p m is the mean 
density of the working gas, T m is the cross-section averaged temper¬ 
ature of the fluid, and r = VT m /T cnt is the ratio of actual tempera¬ 
ture gradient VT m to the standing wave critical gradient 
Tcrit = T m PP\l p m CpU\. 

As shown in Eqs. (1) and (2), the heat power and the produced 
acoustic power are proportional to the temperature gradient factor 
(T - 1). When VT m = VT cr u (i.e. (T - 1) = 0), there is neither heat 
power transferred nor acoustic power produced. For VT m > VT crit 
(i.e. (T - 1)>0), the heat power is transferred away from the 
PAN, and the acoustic power is produced in the regenerator. For 
VT m < VT crit (i.e. (r - 1) < 0), the heat power is transferred to- 


PAN 



Fig. 1. In the standing wave, the regenerator is located on the left side of the PAN. 
The positive direction of x is taken the rightward direction. 


ward the PAN, and the acoustic power is consumed in the regener¬ 
ator. There are two cases for VT m < VT crit . In the case of 
0 < VT m < VT C nt, the regenerator consumes acoustic power to 
pump heat from the cold end to the hot end. In the other case of 
VT m < 0, the regenerator dissipates acoustic power, and the heat 
power is transferred from the hot end to the cold end. In that case, 
the regenerator will not realize the conversion from heat to acous¬ 
tic energy, and there is also no heat pumping from the cold end to 
the hot end. The regenerator located on the right side of the PAN 
can be analyzed in the same way. 

Therefore, some conclusions can be obtained as shown in Fig. 2. 
The required working conditions for the standing wave engine are: 
(1) the hot end of the regenerator is close to the PAN; (2) 
VT m > VTcrit on the left side of the PAN, or VT m < -VT cr n on the 
right side of the PAN. The required working conditions for the 
standing wave cooler are: (1) the ambient end of the regenerator 
should be close to the PAN; (2) 0 < VT m < VT cnf on the left side 
of the PAN, or -VT cr n < VT m < 0 on the right side of the PAN. 


2.2. Traveling wave mode 

The traveling wave thermoacoustic device was originally pro¬ 
posed by Ceperley in 1979 [3]. In the traveling wave heat engine, 
the heat flow direction is opposite to the direction of wave propa¬ 
gation: the wave propagates from the cold end to the hot end, 
while the heat transfers from the hot end to the cold end. The 
regenerator converts heat to acoustic energy when the wave prop¬ 
agates from the cold end to the hot end. In a traveling wave heat 
pump, the wave pumps heat in the opposite direction to its prop¬ 
agation: the wave propagates from the hot end to the cold end, 
while the heat is pumped from the cold end to the hot end. The 
regenerator consumes the acoustic energy to pump heat from the 
cold end when the wave propagates from the hot end to the cold 
end. Thus, the regenerators of the TAE and the TAC in the traveling 
wave should be arranged according to Fig. 3. 


2.3. Traveling-standing wave mode 

We can analyze the thermoacoustic effect of the regenerator ar¬ 
ranged as the four cases in Fig. 4 by combining the thermoacoustic 
conversion performance of traveling wave and standing wave. 


PAN 



(a) The regenerator arrangements of TAE 


PAN 



(b) The regenerator arrangements of TAC 

Fig. 2. In the standing wave, the better arrangements of regenerators. 
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(a) The regenerator of TAE 


Hot end 



Ambient end 


For the TAC, heats are transferred from the hot end to the cold 
end by both the SWC and the TWC. The regenerator cannot pump 
heat from the cold end. 

Case 3. The hot end is close to the PAN, and the TWC propagates 
from the hot end to the cold end. 


Taveling wave 
(b) The regenerator of TAC 



Cold end 
Taveling wave 


For the TAE, when VT m > VT cn t, the SWC converts heat to 
acoustic energy, while the TWC consumes acoustic energy. When 
the thermoacoustic effect induced by the SWC is stronger than that 
by the TWC, the regenerator can produce acoustic energy. 

For the TAC, when 0 < VT m < VTcnt, the regenerator can pump 
heat from the cold end by both the SWC and the TWC, which can 
improve the thermoacoustic gain and efficiency. 

Case 4. The cold end is close to the PAN, and the TWC propagates 
from the hot end to the cold end. 


Fig. 3. In the traveling wave, the better arrangements of regenerators. 

Case 1. The hot end is close to the PAN, and the TWC propagates 
from the cold end to the hot end. 

For the TAE, when VT m > VT cnf , the regenerator realizes the 
conversion from heat to acoustic energy by both the SWC and 
the TWC, which can improve the thermoacoustic gain and 
efficiency. 

For the TAC, when 0 < VT m < VT crit , the SWC pumps heat from 
the cold end to the hot end, while the TWC transfers heat from the 
hot end to the cold end. The direction of the heat transfer induced 
by SWC is opposite to that induced by TWC. When the thermoa¬ 
coustic effect induced by the SWC is stronger than that by the 
TWC, the regenerator can pump heat from the cold end. 

Case 2. The cold end is close to the PAN, and the TWC propagates 
from the cold end to the hot end. 

For the TAE, the SWC consumes acoustic energy, while the TWC 
converts heat to acoustic energy. When the thermoacoustic effect 
induced by the TWC is stronger than that by the SWC, the regener¬ 
ator can produce acoustic energy. 


( a ) PAN 




PAN 




Taveling wave 


Fig. 4. The four cases of the regenerator arrangements. 


For the TAE, both the SWC and the TWC consume acoustic en¬ 
ergy. The regenerator cannot realize the conversion from heat to 
acoustic energy. 

For the TAC, the SWC transports heat from the hot end to the 
cold end, while the TWC pumps heat from the cold end to the 
hot end. When thermoacoustic effect induced by the TWC is stron¬ 
ger than that by the SWC, the regenerator can pump heat from the 
cold end. 

To sum up, the regenerator arranged as Case (4) just consumes 
acoustic power. The regenerators arranged as Cases (1), (2) and (3) 
can produce acoustic power under certain conditions. Especially, 
Case (1) is the best arrangement for the TAE’s regenerator in the 
four cases, because the regenerator of TAE can realize the conver¬ 
sion from heat to acoustic energy by both the SWC and the TWC. 
The regenerator arranged as Case (2) cannot pump heat from the 
cold end, while the regenerators arranged as Cases (1), (3) and 
(4) can pump heat under certain conditions. Especially, Case (3) 
is the best arrangement for the TAC’s regenerator of the four cases, 
because the regenerator can pump heat from the cold end by both 
the SWC and the TWC. 


3. Traveling-standing wave 


Consider wave propagating in the x direction in an ideal gas 
within a channel with constant cross-sectional area A. The hydrau¬ 
lic radius is conventionally defined as [11] 

r h = A/ n. (3) 


In the resonator tube, d v ^:r h , S k < r h . So, in a channel with large 
cross-section, the viscosity due to the interaction between the 
acoustic wave and the solid channel wall can be neglected. With 
the first order of acoustic approximation, the x direction component 
of the momentum equation is 



dU\ 

~dt 


9P± 

dx ’ 



In the ideal acoustic channel, the acoustic field can be described 
by the superposition of a rightward traveling wave with a leftward 
traveling wave. The oscillating pressure can be expressed as 


p 1 = tp + + (1 - T)p_ = p a [ze l{cot kx) + (1 - t)e !( ® t+/o °] 
= p a e icot [ze~ ikx + (1 - T)e ikx ], 



where k is the wave vector, p a is the pressure amplitude, z is defined 
as the traveling wave ratio, and 0.0 ^ z ^ 1.0. z = 0.0 means the pure 
leftward traveling wave, z = 0.5 means the pure standing wave, 
t = 1.0 means the pure rightward traveling wave, 0.0 < z < 0.5 
means the mixed acoustic field of the SWC with the ratio of 2 t 
and the leftward TWC with the ratio of (I-t)-t, and 
0.5 < t < 1.0 means the mixed acoustic field of the SWC with the ra- 
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tio of 2(1 — t) and the rightward TWC with the ratio of 
t - (1 - ^.Substituting Eq. (5) into Eq. (4) 


Ul 


1 

Pm 


$! dt = R^ [Te - to _ 

ax p m c 


(1 - T)e Ikx ], 



where c is the sound speed. 

According to Eqs. (5) and ( 6 ), the specific acoustic impedance 
can be yield, 


„ Pi Te~ lkx + (1 - T)e lkx 

Ui Te~ !/<x - (1 - T)e lkx 



The phase of Z represents the leading phase cp z of p i to U\. Fig. 5 
shows the results of the calculation for (p z versus the normalized po¬ 
sition kx/jL for different t. When t = 0.0, it is the pure leftward trav¬ 
eling wave, and (p z = ± 180°. When t = 0.5, it is the pure standing 
wave, and (p z = -90° in the region of - 7 i /2 <kx< 0 and cp z = 90° in 
the region of 0 </<x< 7 i/ 2 . When t = 1.0, it is the pure rightward 
traveling wave, and cp z = 0°. When 0.0 < t < 0.5, it is the mixed 
acoustic field of the SWC and the leftward TWC, and -180° < 
(p z < -90° in the region of —tt/2 <kx< 0 and 90° < cp z < 180° in the 
region of 0 < kx < tt/ 2. When 0.5 < t < 1.0, it is the mixed acoustic 
field of the SWC and the rightward TWC, and -90° < cp z < 0° in the 
region of —zr/2 < kx < 0 and 0° < cp z < 90° in the region of 0 < kx < n\2. 

According to the analysis and Fig. 4, it is obtained that -90° < 
(p z < 0° for Case (1), 0°<^ z <90° for Case (2), -180° < (p z < -90° 
for Case (3), and 90° < (p z < 180° for Case (4). 


4. Thermoacoustic formula analysis 

This section discusses the thermoacoustic effect of one point of 
the regenerator instead of the whole regenerator. Although, the 
acoustic field cannot keep invariant along the regenerator, the 
acoustic field in the regenerator can be regarded as a set of discrete 
points with variant acoustic characteristics. We will locate the 
regenerator in a certain region of acoustic field, which has better 
characteristics for thermoacoustic conversion. 


4A. Acoustic power 


The average acoustic power dE 2 produced in the length dx of a 
regenerator is expressed as [ 11 ], 


dE 2 _ 1 cop m lm[-/, 

dx ~ 2 A 


|1 -fv\ 


2 — 2 1 

2 y p m 2 


i dr 


X 


m 


T m dx 


Re 


fk-fv 


.(1 -W(i-ff) 


Pi l/, 


( 8 ) 


N 


180- 

120 - 

60- 

0 - 

-60- 

- 120 - 

-180 - 
-0 


X A°°Oooooo0OOOOOO 0 ^A 

A a a A 

a Maaaaaa aA^ 

©©©©©©©©©©©©©©©©©©©© 


A. .A 

A A AaaaaaAA A 

©©©©©©©©©©©©©©©©©©©© 

aAAAAAAZVVv\^ 

a A 

/ooo oOOQOOOOo Oo 0 A ' a 

Ai 08 nnnnnnnnnnnnnnn 80 C 


a aaaaaaaa Aa 


Mi 


A 

„ ▲ 

:i« 


.6 


-0.4 


- 0.2 


0.0 

kx/71 


0.2 


0.4 


— ■ — x=0.0 

-•—i=0.2 

-a— x= o .4 

— © — t= 0.5 
A x=0.6 
o x=0.8 
n x=l .0 


0.6 


Fig. 5. cp z versus kx/n with different t: 0.0, 0.2, 0.4, 0.5, 0.6, 0.8, 1.0. 


where cr is the Prandtl number, 7 is the ratio of isobaric specific heat 
to isochoric specific heat, p m is the mean pressure of the fluid, Ui is 
the volume flow rate, and / v , f k are the cross-sectional averaged 
Rott’s functions which depend on the specific channel geometry un¬ 
der consideration. 

The functions of f v and f k are known for many geometries. As 
shown in Fig. 5.14 of Ref. [11], the solutions for most of these 
geometries are actually very similar. Therefore, this paper only 
takes the parallel-plate type of regenerator as a case to investigate, 
and other types of regenerator can be analyzed in the similar way. 
y = 0 is defined at the centerline of the spacing between two adja¬ 
cent parallel plates, 


tan h[(l +i)y Q /d v \ 
(1 + OTo/<5 v 
tan h[(l +i)yjd k ] 
(i + i)yj dk 


( 9 ) 

( 10 ) 


where 2 y 0 is the thickness of each spacing, and r h = y 0 . 

The relative hydraulic radius is defined as £ = r h lS kl the ratio of 
the hydraulic radius r h to the thermal penetration depth S k . The 
functions of f v and f k can be written as follows, 



fk 


tan h[(l +i)Cy / V^^ 

(1 + i)C\/T ~Jg 

tan h[(l +i)C] 

(1 + Of ■ 


( 11 ) 

( 12 ) 


The first two terms in Eq. ( 8 ) are always negative, which con¬ 
sume the acoustic power, and are independent of temperature gra¬ 
dient along x. The first term describes the viscous dissipation of 
sound wave. It is proportional to |Ffi | 2 and independent of thermal 
conductivity. The second term describes the thermal-relaxation 
dissipation. It is proportional to |pi | 2 and is independent of viscos¬ 
ity. The third term is called as the source/sink term, because it can 
either produce or consume acoustic power. It exists only if dT m / 
dx 7 ^ 0. For thermoacoustic engine, it is positive because the ther¬ 
mal energy is converted to the acoustic energy. It can be written as 



1 ( 1 dT m \ 

2 \T m dx) 


|Pil|Ui||e g |cos(<p e 




where 



fk -fv 

(1 —/»)(! - ff) 


e g \e 


We 



and (p e is the phase of e g . 

As shown in Eqs. ( 8 ) and (13), the maximum power gain dE 2 /dx 
of the TAE appears at (p z = (p e when dE m \dx > 0. Fig. 6 shows (p e ver¬ 
sus C, and -64° ^ q> e < 0 °. So, the optimal leading phase appears in 
the region of -64° ^ (p z = (p e (£) ^ 0°, corresponding to the Case (1) 
of Fig. 4. 

As shown in Fig. 6 , with the increase of £, \e g \ decreases. For the 
idea regenerator without viscous dissipation, when £ -> 0 , (p e -> 0 ° 
and \e g \ reaches the maximum \e g \ 1.0. So, dE 2 /dx reaches the 
maximum at the traveling wave phase (p z = (p e = 0°. However, for 
the practical regenerator, the viscous dissipation cannot be ne¬ 
glected, and sometime is the main dissipation. So the moderately 
small hydraulic radius is adopted in the practical traveling wave 
device. It leads (p e to deviate from 90° by thermo-relaxation, which 
indicates that the optimal leading phase deviates from the travel¬ 
ing wave phase, (p z = (p e ¥^ 0°. So the optimal q> z is not the traveling 
wave phase but traveling-standing wave phase with a negative va¬ 
lue — 64° ^ (p z = (PeiO ^ 0°. 

In order to yield confidence in the analysis, we now give some 
numerical results about acoustic power amplification comparing 
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Fig. 6. \e g \ and cp e versus £. 



with the experimental data of Biwa et al. [4] who measured the 
acoustic power amplification under different cp z . 

Acoustic power amplification can be expressed as, 


Eout _ Ein + A E _A E _ 

Ein - E in ~ + 0.5|p 1 ||Ui | cos (cp z ) ’ 


Substituting Eq. (8) into Eq. (15) 


Eout ^ 271/ 

E in 2cos(c p z ) 


-j^xpL-(y-l)Im[-/, ( ]x|Z N 
+ ^iifl e gl cos((p e -(p z ) 



According to the experiment of Ref. [4] and Eq. (16), the acoustic 
power amplifications versus cp z are calculated when £ = 0.36 and 
1.87, as shown in Fig. 7. The theoretical curves have the similar 
shapes as corresponding experimental curves, showing a general 
agreement between the theoretical and experimental results. Thus, 
we can explain the experimental results with the viewpoint dis¬ 
cussed above. 

As shown in Fig. 7a, when A T m > 0, E 0Ut IE in shows a maximum 
reaching 1.7 at cp z « -42°.This corresponds to the situation where 
the hot end of the regenerator is located close to the velocity node. 
Biwa et al. [4] explained this phenomenon as follows: “Since the 
kinematic viscosity /i of the gas increases with increasing temper¬ 
ature, viscous losses should become relatively larger at the hot end 
than that at the cold end, even if the velocity amplitude is the same 
at both ends of the regenerator. Thus, <D m = - cp z yielding minimum 
viscous losses should be positive so that the velocity node comes 
close to the hot end. Therefore, we consider that the gain shows 
a maximum at <D m = - cp z > 0° instead of <D m = - (pcp z = 0°.” How¬ 
ever, calculations in this paper use the average values of regenera¬ 
tor parameters, and neglect the changing of /a with the 
temperature. The maximum acoustic power amplification also ap¬ 
pears at <D m = - cp z > 0°. That reason is that cp e deviates from 0° to a 
certain value between 0° and -90° because of £ = 0.36. So the opti¬ 
mal value of -90° < (p z < 0° is chosen to match the deviation of cp e . 
When A T m > 0 and cp z < 0°, this corresponds to the arrangement as 
Case (1) in Fig. 4, which is the best arrangement for the TAE. The 
regenerator realizes the conversion from heat to acoustic energy 
by both the SWC and the TWC. As shown in Fig. 7a, when 
A T m < 0, the maximum of E 0Ut IE in appears in the region of cp z > 0°. 
For this phenomenon, Biwa et al. [4] gave the similar explanation 
for the phenomenon when A T m < 0: “<D m = -cp z yielding minimum 
viscous losses should be negative so that the velocity node comes 
close to the hot end.” However, calculations in this paper use the 
average values of the regenerator’s parameters, and neglect the 
changing of /a with the temperature. The maximum acoustic power 




Fig. 7. E 0Ut IE in versus cp z for a variety of A T m and £: (a) £ = 0.36, (b) £ = 1.87. 


amplification also appears at <D m = - 4> z < 0°. That reason is that 
cpcp e deviates from 0° to a certain value between 0° and -90° be¬ 
cause of £ = 0.36. Thus the value of cos (0 e - cpcp z ) = cos cp e cos cp z + - 
sin cp e sin (p z in the region of 90° > cp z > 0° is smaller than that in the 
region of 0° > 0 Z > - 90°. This indicates that the acoustic power is 
smaller in the region of 90° > cp z > 0°. 

As shown in Fig. 7b, the curve of E 0 Ut IE in versus cp z when £ = 1.87 
is clearly different from that when £ = 0.36. When A T m > 0, E 0 Ut IE in 
monotonically decreases with the increase of cp(p z and it exceeds 
T h IT c = 1.9 when 0 Z < -60°. When £ = 1.87, the thermal-relaxation 
dissipation becomes the main loss, and the acoustic power gain 
is mainly dependent on the thermoacoustic conversion of the 
SWC. In the standing wave engine, E in 0.0 and acoustic energy 
can be generated within the regenerator. This represents that 
Eout IE in -> oo when the energy conversion is mainly dependent on 
the SWC. Thus, the results clearly demonstrate that amplification 
of E 0 Ut IE in exceeding T H IT C is possible through the SWC energy con¬ 
version in the regenerator with £ = 1.87 [4]. 

Fig. 7 can also be explained through the analysis of Section 2 in 
this paper. When A T m > 0 and qxp z < 0°, the regenerator arrange¬ 
ment is the same as Case (1) of Fig. 4: The regenerator realizes 
the conversion from heat to acoustic energy by both the SWC 
and the TWC. Thus, E 0 Ut IE in > 1.0 in the region of cp z < 0°, as shown 
in Fig. 7. When A T m > 0 and cp z > 0°, the regenerator arrangement is 
the same as Case (2) of Fig. 4: The SWC consumes acoustic energy, 
while the TWC converts heat to acoustic energy. When the ther¬ 
moacoustic effect induced by the TWC is stronger than that by 
the SWC, the regenerator can produce acoustic energy. Thus, F out / 
E in >1-0 appears in the region of cp z > 0°, as shown in Fig. 7. When 
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the thermoacoustic effect induced by the SWC is stronger than that 
by the TWC, the regenerator contributes to the damping of acoustic 
energy. Thus, E 0Ut IE in <1.0 also appears in the region of q> z > 0°, as 
shown in Fig. 7. When A T m < 0 and cp z > 0°, the regenerator arrange¬ 
ment is the same as Case (3) of Fig. 4: the SWC converts heat to 
acoustic energy, while the TWC consumes acoustic energy. When 
the thermoacoustic effect induced by the SWC is stronger than that 
by the TWC, the regenerator can produce acoustic energy. Thus, 
E 0 utlEin> 1.0 appears in the region of cp z > 0°, as shown in Fig. 7. 
When the thermoacoustic effect induced by the TWC is stronger 
than that by the SWC, the regenerator contributes to the damping 
of acoustic energy. Thus, E 0Ut IE in <1.0 also appears in the region of 
cp z > 0°, as shown in Fig. 7. When A T m < 0 and cp z < 0°, the regener¬ 
ator arrangement is the same as Case (4) of Fig. 4: both the SWC 
and the TWC consume acoustic energy. Thus, E 0Ut IE in < 1.0 in the 
region of cp z < 0°, as shown in Fig. 7. 

According to the analysis above, the optimum value of cpcp z is 
not only dependent on the changing of /i with the temperature, 
but also mainly dependent on the phase deviation because of the 
thermo-relaxation with the finite hydraulic radius. The optimum 
value of (p z must match the phase deviation due to the thermo¬ 
relaxation. It is found that the optimum traveling-standing wave 
phase cp z can also make both the TWC and the SWC contribute to 
the thermoacoustic conversion. 

4.2. Heat flow 



Fig. 8. cp q versus £■ 
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In the regenerator model, ordinary thermal conduction in x 
direction is neglectable. The average heat flow generated along 
regenerator due to hydrodynamics transport is [1,12] 
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When dT m \dx > 0, the first term in Eq. (17) is negative, which trans¬ 
ports undesired heat from hot end to cold end, and is independent 
of cp z . The first term describes the acoustically induced heat convec¬ 
tion. The second term in Eq. (17) can either pump heat from cold to 
hot or transport heat from hot to cold. For the TAC, it is positive 
when dT m ldx > 0, because the acoustic wave pumps heat from cold 
to hot. It can be written as, 

Q.g = - 0 IPi IIIA |) |<j g | cos (cp q + cp z ), (18) 

where 



fk -fv 

(1 + <7)(1 -/„) 




and (p q is the phase of q g . 

As shown in Eq. (18), the maximum cooling power Q of the TAC 
appears at q> z = 180°- cp q when dT m \dx > 0. Fig. 8 shows cp q versus £, 
and -84° < cp q < 0. So, the optimal leading phase appears in the re¬ 
gion of -180° < cp z = 180 °-cp q < -96°, corresponding to the Case (3) 
of Fig. 4. 


4.3. Second law efficiency of the TAE 


The second law efficiency r} 2 is the first law efficiency y \i divided 
by the Carnot efficiency rj c , 


=Pl/Pr 


According to Eqs. (8), (17), and (20), it can be yielded, 


where \j/ is the operation factor of the regenerator, defined as 


c 1 dT m 

co T m dx ’ 


( 22 ) 


q 2 has been calculated for a range of cp z and £ according to Eq. 
(21) on a minicomputer having a FORTRAN complier capable of 
evaluating complex equations. In the calculations, we use the 
values of the thermoacoustic Stirling heat engine described in 
the Ref. [13], such as |Z N | = 15, if = 22 and helium as working 
gas. 

Fig. 9 shows the relation between r/ 2 and cp z with different £. 
With the increase of £, the maximal r\ 2 of each £ increases at first, 
reaches the maximal value 0.786 at the point of (£ = 0.21, 
cp z = -25.2°), and then decreases. In all cases of £, rj 2 reaches the 
maximum in the region of -90° < cp z < 0, which is relative to Case 
(1) in Fig. 4. When £ -► 0, the optimal cp z is close to 0°. However, 
rj 2 0, because viscous losses should become infinite when 



( 20 ) 


Fig. 9. rj 2 versus cp z for a variety of (. The fine solid line connects the optimal points 
of different £. 































456 


H. Kang et al.JCryogenics 50 (2010) 450-458 




Fig. 10. (a) The maximum r/ 2 with different i J/. (b) The optimal £ and the optimal cp z 
for the maximum r\ 2 with different iff. 




Fig. 11. (a) The maximum rj 2 with different \Z N \. (b) The optimal £ and the optimal 
cp z for the maximum r\ 2 with different \Z N \. 


£ -► 0. With the increase of £, the optimal cp z deviates from 0° to a 
certain value between 0° and -90°. 

Furthermore, the maximum r\ 2 , the optimal £ and the optimal cp z 
with different i jj are summarized in Fig. 10. With the increase of i J/, 
the maximum r\ 2 increases sharply at first, and then increases 
slowly and gradually reaches to r\ 2 = 0.85. As shown in Fig. 10b, 
as the optimal £ deviates from 0.0, the optimal cp z deviates form 
0° because of the thermo-relaxation. Although the range of the 
optimal (p z is large, all of the optimal cp z is located in the range of 
-90° to 0°. 

On the other hand, the maximum r/ 2} the optimal £ and the opti¬ 
mal cp z with different \Z N \ are summarized in Fig. 11. Similarly to 
the phenomenon of i j/ } with the increase of \Z N \ } the maximum r\ 2 
increase, and then reaches to r\ 2 = 0.94. The optimal phase differ¬ 
ence cp z depends on the optimal £ because the phase difference 
needs to balance the phase deviation generated by thermo-relaxa¬ 
tion. Similarly, all of the optimal cp z is located in the range of -90° 
to 0°. 

As a result, the TAE’s regenerator arranged as Case (1) is the 
best. 

4.4. Second law coefficient-of-performance of the TAC 

The second law coefficient of performance COP 2 is the first law 
coefficient of performance COP! divided by the Carnot coefficient 
of performance COP c 

COP 2 = COP!/COP c . (23) 


According to Eqs (8), (17), and (23), 
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(24) 



Fig. 12. COP 2 versus cp z for a variety of £. The fine solid line connects the optimal 
points of different £. 
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Fig. 13. (a) The maximum COP 2 with different i j/. (b) The optimal £ and the optimal 
cp z for the maximum COP 2 with different \Jj. 




Fig. 14. (a) The maximum COP 2 with different \Z N \. (b) The optimal £ and the 
optimal cp z for the maximum COP 2 with different \Z N \. 


COP 2 has been calculated for a range of cp z and £ values according to 
Eq. (24) on a minicomputer having a FORTRAN complier capable of 
evaluating complex equations. In the calculations, we let \Z N \ = 15, 
i/j = 3.6 and use helium as working gas according to the thermoa¬ 
coustic Stirling cooler described in Ref. [14]. 

Fig. 12 shows the relation between COP 2 and cp z with different £. 
With increase of £, the maximal COP 2 of each £ increases at first, 
reaches the maximal value 0.85 at the point of (£ = 0.24, 
(p z = -149.4°), and then decreases. In all cases of £, COP 2 reaches 
maximum in the region of -180° < cp z < -90°, which is relative to 
the Case (3) in Fig. 4. When £ -> 0, the optimal cp z is close to 
-180°. However, COP 2 -► 0, because viscous losses should become 
infinite when £ -► 0. With the increase of £, the optimal cp z deviates 
from -180° to a certain value between -180° and -90°. 

Furthermore, the maximum COP 2 , the optimal £ and the optimal 
(p z with different i/^ are summarized in Fig. 13. With the increase of 
i j/, the maximum COP 2 increase sharply at first, and then increase 
slowly and gradually reaches to COP 2 = 0.85. As shown in 
Fig. 13b, as the optimal £ increases from 0.0, the optimal cp z in¬ 
creases form -180° because of the thermo-relaxation. Although 
the range of the optimal cp z is large, all of the optimal cp z is located 
in the range of -180° ~ -90°. 

On the other hand, the maximum COP 2 , the optimal £ and the 
optimal cp z with different \Z N \ are summarized in Fig. 14. Similarly 
to the phenomenon of i j/, with the increase of \Z N \, the maximum 
COP 2 increase, and then reaches to COP 2 = 0.73. The optimal phase 
difference cp z depends on the optimal £ because the phase differ¬ 
ence needs to balance the phase deviation generated by thermo¬ 
relaxation. Similarly, all of the optimal cp z is located in the range 
of -180° to -90°. 


As a result, the TAC’s regenerator arranged as Case (3) is the 
best. 

5. Conclusions 

Actually, the acoustic field in practical thermoacoustic devices 
is neither a pure traveling wave nor a pure standing wave. The 
thermoacoustic effect is the hybrid effect of traveling wave compo¬ 
nent and standing wave component. Thus, the thermoacoustic ef¬ 
fect of the traveling-standing wave is investigate in this paper. 
Firstly, the thermoacoustic conversion properties of the travel¬ 
ing-standing wave are analyzed qualitatively by combining the 
thermoacoustic conversion properties of the traveling wave with 
those of the standing wave. Then, according to the basic thermoa¬ 
coustic formulas, this paper analyzes the optimum phase differ¬ 
ence which can make both the traveling wave and the standing 
wave contribute to the thermoacoustic conversions. The results 
are consistent with the qualitative analysis: to gain a higher acous¬ 
tic power and efficiency, the hot end of the regenerator in a ther¬ 
moacoustic engine should be close to the pressure antinode, and 
the traveling wave component should propagate from the ambient 
end to the hot end. To gain a higher cooling power and coefficient 
of performance, the ambient end of the regenerator in a thermoa¬ 
coustic cooler should be close to the pressure antinode, and the 
traveling wave component should propagate from the ambient 
end to the cold end. 

Furthermore, this paper analyzes the influence of the parame¬ 
ters of the acoustic field and the regenerator’s structure on the 
thermoacoustic conversion, and finds out the optimum condition 
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for the thermoacoustic conversion. The optimal phase difference 
between pressure and velocity depend on the optimal hydraulic ra¬ 
dius because the phase difference needs to balance the phase devi¬ 
ation generated by thermo-relaxation. 

Additionally, our theoretical results also show a good agree¬ 
ment with the experimental data [4], which validates the analysis 
in this paper. Besides, the analysis in this paper further provides a 
more intensive understanding of these experimental results. For 
improving thermoacoustic conversion, the optimum cp z is not only 
dependent on changing of the kinematic viscosity with the temper¬ 
ature, but also mainly dependent on the phase deviation because of 
the thermo-relaxation with the finite hydraulic radius. 
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